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Abstract 

This Proceeding presents the method that allows us to get the Gyro-Kinetic Approximation of the 
Dynamical System satisfied by the trajectory of a particle submitted to a Strong Magnetic Field. The 
goal of the method is to build a change of coordinates in order to make the dynamic of two components of 
the trajectory to disappear. This change of coordinates is based on a Darboux mathematical Algorithm 
and on a Lie Transform 



1 Introduction 

Scientific framework - At the end of the 70', Littlejohn [TTJ [TU shed new 
hght on what is cahed the Guiding Center Approximation. His approach was based 
on mathematical theories - Hamiltonian Mechanics, Differential Geometry, Symplectic 
Geometry - in order to clarify what has been done for years in the domain (see Kruskal 
[TOJ, Gardner [4J, Northrop [14J, Northrop & Rome [15j). His papers claim that it has 
been for him an enormous effort to reach this goal, since he had to incorporate into a 
physical affordable theory high level mathematics. Sure, this theory is a nice success. 
It has been being widely used by physicists to deduce related models {Finite Larmor 
Radius Approximation, Drift Kinetic Model, Quasi-Neutral Gyro-kinetic Model, etc.^ 
see for instance Koseleff [9], Brizard [T], Dubin et al. [2], Frieman & Chen [3j, Hahm 
[7], Hahm, Lee & Brizard [8j, Parra & Catto [16) |17l [18]) making up the Gyro-Kinetic 
Approximation Theory^ which is the basis of all kinetic codes used to simulate Plasma 
Turbulence emergence and evolution in Tokamak (see for instance Grandgirard et al. 

EE]). 

The incorporation of mathematical concepts into a physical theory has been done so 
nicely that the resulting Gyro-Kinetic Approximation Theory is now very difficult for 
mathematicians. 
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This talk (and proceeding) is the first step into the rewriting of the Gyro-Kinetic 
Approximation Theory into a mathematicahy affordable theory. It is only a summa- 
rize, to understand, with a mathematical slant, several aspects of what is done in the 
references just cited. 

Charge particles submitted to Strong Magnetic Field - The context of 
the Gyro-kinetic Approximation is Tokamak Physics. An artist vision of Iter, which 
is a Tokamak, is given in Figure [T] The vessel of a Tokamak is the interior of a torus 
with a vertical axis of symmetry. Along the torus, electromagnets can generate a large 
magnetic field. 




Figure 1: Artist vision of Iter. 



We begin by considering a charged particle within the vessel of a Tokamak in service. In 
Usual Coordinates (x, v) = (xi, X2, X3, 1'l, '^3)7 where x stands for the position vari- 
able and V for the velocity variable, the position and the velocity at time t (X(t; x, v, s), 
V(t; X, V, s)) of the particle which is in x with velocity v at time s is the solution to 

V, (1) 

^(E(X)+VxB(X)). (2) 

m 

This dynamical system means nothing but that time derivative of the position is the 
velocity and that time derivative of the velocity, which is the acceleration, in linked with 
Lorentz force by Newton's law. In this dynamical system magnetic field B is composed 
of a strong applied piece, of a strong self induced piece and of a self induced pertur- 
bations. The strong applied piece is generated by the electromagnets and wind around 
the axis of symmetry of the tokamak. Charged particles winding, with a relatively 
large velocity, around the axis of symmetry within the vessel generate the indispensable 
vertical component of the large magnetic field. This is what is called here the strong 
self induced piece. The part of the self induced magnetic field which is not vertical is 
called here the self induced perturbations and is, from now, forgotten. Electric field E 
is induced by the particles of the vessel. 

Without loss of generality, we can consider that E is the opposite of the gradient of 
an electric potential, and B the curl of vector potential, i.e.: 

E = -V<l>, B = V X A. (3) 

Helicoidal trajectories - Larmor Radius - It is weh known that the trajec- 
tory of a charged particle submitted to a magnetic field is a helix. The helix axis is 
the magnetic field direction and its radius, which is called in the context of Tokamak 



ax 

dt 
d\ 

~dt 
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Plasma the Larmor Radius, equals the norm of the projection of the velocity on the 
plan orthogonal to the magnetic field divided by the particle's mass time the norm of 
the magnetic field. As illustrated in left picture of Figure [2j because of the mass ratio 




Figure 2: Left: Helicoidal trajectories of ions and electrons submitted to a magnetic field. 
(Source: S. Jardin's Lectures at Cemracs'lO). Right: Screen effect and Debye Length in a 
plasma. 



between ions and electrons, the order of magnitude of electron Larmor Radius is much 
smaller than that of ions. In Tokamak, those order of magnitude are ^ 5 • 10~^m for 
electrons and ^ 10~^m for ions. Those two lengths are two important scales of Tokamak 
Plasma Physics. 

To simplify the purpose, from now we will forget the scale of electron Larmor Radius 
and only take ions under consideration. 

Debye Length - We now explain an important phenomenon of Plasma Physics 
which is the Screen Effect, with its related length which is the Debye Length. 

It is well known that a charged particle generates around it an electric potential 
which depends on the inverse of the distance to it, say ^ 1/r. Since this function is 
the Green Kernel of the opposite to the Laplace operator, when following the Mean 
Field routine to deduce the model taking under consideration interactions of many 
particles, it is gotten that the electric field that fills the space is the solution to —A^ = 
Cst Charge Density. Because, of the long range of ~ (l/r)-potential, the regularization 
properties of minus Laplacian are such that space variations of the electric potential (and 
consequently of the electric field) occure at a large scale. 

Now, if the density of particles is relatively large, as it happens in Plasma, and 
if we consider a test particle (which is drawn in blue in right picture of Figure |2|, 
the other particles (in the picture, orange particles are positively charged and green 
particles negatively charged) which are beyond a given length A, which is the Debye 
Length (the black circle in the picture is censed being of radius A), may be gathered 
into subsets (the two grey ellipses in the picture are such subsets) which are such that 
the resulting action of pairs of those subsets on the test particle is negligible. This has 
the following consequence: it can be considered that the potential ~ 1/r generated by 
a particle may be approximated by a ~ e~^^^/r potential. Then, following the Mean 
Field routine with this new assumption leads an electric field that fills the space which 
is the solution to —V^ = Cst Charge Density. Operator —V is the pseudo-differential 
operator which Green Kernel is e~^/^/r whose regularization properties are much less 
than that of minus Laplacian. Then the regularity of the electric field is not as high as 
we could think before this formal analysis. In particular the electric field can be prone 
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to variations over lengths ranging from several Debye Lengths to several hundreds of 
Debye Lengths. Those variations are maybe not extremely large, but as no principle 
fends them off, they are stable. 

What Ions see - In a very simplified slant. Electric field in a plasma is a bed with 
variations which are not so large but not negligible, and, more than anything, stable 
at scales ranging from several Debye Lengths to several hundreds of Debye Lengths. In 




Figure 3: Symbolic picture of oscillating trajectory of ion in a electric field with variations 
with typical size comparable with trajectory oscillation. 

Figure [3] we symbolically represent this bed as a charcoal grey surface. (It is symbolic 
because a surface is drawn while, in real world. Electric Field fills a 3D space.) This 
surface shows variations at two scales which are cenced to be about 100 times the Debye 
Lengths and 1.000 times the Debye Lengths. Those variations are not of large amplitude 
but are stable. Now, ion motion, symbolically represented in red in the figure, shows 
an oscillation of amplitude which is about 100 times the Debye Length (ion Larmor 
Radius ~ 10~^m, Debye Length : ~ 10~^m). Hence, because of the stability of bed 
variations, its resulting effect can be important on ion trajectory after several periods 
of oscillation. To take into account this resulting effect, the concept of gyro-average was 
introduced. 

Dimensionless Dynamical System - Here, we do not present the scaling rou- 
tine, we only mention that two small parameters may be introduced. The first one 
is 

Ion Larmor Radius 10~^m o 

e ~ 10 (4) 

Tokamak size 10m 

The second one is 77 the Debye Length linked variations of the Electric Field and Po- 
tential. It is defined as 

Characteristic scale of variations of the Electric Field . . 

Tj ~ . (5) 

Tokamak size 
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If we follow a scaling routine, we will get that the dimensionless Electric Potential and 
Field write 

$(x) = $o(x)+7?^i(-) and E(x) =Eo(x)+Ei(-), (6) 

T] T] 

with 

E(x) = -(V$)(x), (7) 
and that the dimensionless Magnetic Field write 

^ = (V X ^)(x). (8) 

Then, the dimensionless trajectory, is solution to the following dynamical system: 

^=Eo(X) + Ei(^)+VxM. (10) 

We guess that anybody agree with the fact that this system is a good model to describe 
motion of charged particle within Tokamak. 

Gyro-kinetic Model - Yet, there is the Gyro-Kinetic Model. It claims that a 
trajectory (r, ^0, , J, F) in a coordinates system (r, ?/^, (/:), j, 7) we will discuss 

on later on, is solution to: 

dt ^ dt q{r)R r ' 

where 




Figure 4: Position coordinates in use in a Tokamak. (Source: V. Grandgirard's Lectures at 
CemracsTO). 

(12) 
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and 



We now explain how to read equations (11) - (13). A Tokamak is the interior of a 



torus with a vertical symmetry axis. (A Tokamak and its associated position coordinate 
system {r^ip^cp) is given in Figure [4j) Hence, it is thought as a small disc carried by 
a large circle of radius Rq. Then, for a given point of the Tokamak, r is (almost) its 
distance from the center of the small disk, is (almost) the angle that the point makes 
in the small disk and (f is (almost) the angle that it makes on the large circle. 
In what concerns the velocity variable system (^^;||, j, 7), if a moving point is in (r, t/j^ cp)^ 
the It; 1 1 -component of its velocity stands (almost) for the projection of its velocity on the 
direction the Magnetic Field has in (r, ip^ (p) (which is not far from being tangent to the 
large circle but which has a small component tangent to the small circle), j is (almost) 
half the square of the norm of its velocity on the plan orthogonal to the Magnetic Field 
direction, and 7 is (almost) the angle that this projection makes in this plan. 

Above, precision almosV^ is essential, and constitutes the core of the Gyro-Kinetic 
Approximation which consists, as we will see, in building infinitesimal changes of coor- 



dinates to set out a dynamical system shaped like (11). 



We now comment system (11), forgetting that the variables are only almost what 



they are supposed to be. The first equation in (11) says that the r-component of a 
particle trajectory varies with a velocity which is the projection on the direction of the 
vector joining the center of the small disk and the point of the Electric Drift velocity ED 



and the Magnetic Curvature Drift velocity MCD^ both defined in (12). The Electric 
Drift velocity ED involves the cross product of the Magnetic Field and what is called 
the gyro-average of the Electric Field, which is indicated by symbol ( ). The Magnetic 
Curvature Drift velocity MCD involves the cross product of the Magnetic Field with 
the gradient of the Magnetic Field norm. 

In the second equation of ( pT| , q(r) is called the quality factor. Essentially, it is the 
number of revolutions in the small disk a particle makes while making one revolution 



along the large circle. Then, term W\\/(q(T)R)^ where R is defined in (12), generates 



the winding of the particle trajectory in the small disk (this is illustrated by the red 



curves drawn on the left of the figure). The second term of the second equation of (11 ) 
is linked with the action of the Electric Drift velocity and Magnetic Curvature Drift 
velocity. 

The third equation involves the contribution W\\/R of the -component of the particle 
velocity to the variation of the (/^-component of the particle position. 
The last equation describes the variation of it; y -component of the particle velocity. It 
involves, among other terms, the projection of the gyro-average of the Electric Field on 
the direction of the Magnetic Field and the variation of the Magnetic Field norm in the 



direction of the Magnetic Field, materialized by operator V\\ defined by (13). 



In system (11), their is neither equation for J (the j-component of the particle 
velocity) nor for F (the 7-component of the particle velocity). This may be considered 
as a strange fact. The reason why this happens is that the dynamics of r,'0, (^, VFy, 
the r-, ?/^-, (^-, It; 1 1 -components of the particle trajectory and velocity do not depend on 
the 7-component of the particle velocity and that the its j-component is invariant. This 



is summarized as: 



dJ dV ^ ,,,, 

-— = 0, -— = Does not matter. (14) 

dt ' dt 
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What we explain and what we do not explain - In this document we give 
the key ideas that bring the deduction of the Gyro-Kinetic Dynamical System ( 12 ) from 
the original one ^ - (fTol). 



Among thing that we do not tackle, there is the question of the specific position 
coordinates {r^tjj^cp). We will remain with a generic position coordinate system. We 
will not talk about Gyro-average either. The question of the coupling of Gyro-Kinetic 
Approximation with computation of Electric and Magnetic Fields, using Quasi-Neutral 
Poisson Equation or Maxwell Equation is also not addressed here. 

To end this paragraph we mention that the method only works in case when 77 is 
larger than e, or in other words, when 77 = 6:^"'^ for /^>0. This is a important limitation 
in the theory. 



2 Methode summarize 



Key result - The method allowing us to go from system d9]) - (10) to a system of 



the kind of ( |ll| ), i.e. to make disappear two equations, has nothing to do with magic. 

(ri, r2, ra, r4, rs, re), a Hamiltonian 



It is based on the following result: 

Theorem 2.1. If, in a given coordinate system r 
Dynamical System writes: 



( 



~dt 



P(R)VriJ(R), P(R) 



M 



0-- 
V 



\ 



-1 



(15) 



(for a trajectory R = i?2, i?3, -R4, R5, Rq)^ )) with a Hamiltonian Function that does 
not depend on the last variable, i.e. 

dH 



drt 



= 0, 



(16) 



0, 



(17) 



Then, suhmatrix Ai does not depend on the two last variables, i.e. 

dM ^ BM 
—— = and —— 

and consequently, time- evolution of the four first components i^2, ^3, R4 is indepen- 
dent of the last component Rq; and, the penultimate component of the trajectory in 
not time- evolving, i.e. 

dRs 



dt 



= 0. 



(18) 



The proof of property (17) is relatively straightforward by a calculation. Equation 



( 18 ) is a direct consequence of ( 16 ) that yields a sixth component of Vri^ (R) with worth 
and consequently, because of the form of V(R) a fifth line of Dynmical System (15) 
which is exactly (18). 



Essentially, the Gyro-Kinetic version of Dynamical System (l9|) - (10) exhibit es- 



sentially a property like (18) and the property that the four first components of the 



trajectory are independent of the last component. 

As a matter of fact, building the Gyro-Kinetic Approximation consists in building a 
change of coordinates that bring a system of coordinates in which 7^(R) has the same 
form as in (15) and in which (16) is true. 
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Usual Coordinates 




(x, 


v) 








Cylindrical Coordinates 


dt 
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4: Darboux Method | 






Darboux Almost 


: Hamiltonian? 




Canonical Coordinates 








Canonical Coordinates 

(q, p) 




5: Lie Method], 


if^- = /i,-(q,p): 
dQ 
dt 






Lie Coordinates 


aP 

'dt ~ 









Usual Coordinates 
(x,v) 



if^-(x,v),7>f(x,v) s.t 

at 

dV 



1: Hamiltonian? 




Canonical Coordinates 

(q, p) 

s.t: 
^dQ\ 

dt 



dP_ 

~dt J 



Cylindrical Coordinates 

{■x,v\\,v_i,9) 
H^{^, v\\,v±,e), (x, v\\,v^,e) 



4: Darboux Method | 





\ 







Darboux Almost 
Canonical Coordinates 



5: Lie Method | 



Lie Coordinates 
(z,^||, j,7) 

5^^(z, , j) , (z, w\\ , j, 7) 



Figure 5: The method is made of 5 steps. 1: Check that the Dynamical System ([9])- (10) is 
well Hamiltonian. 2: Write the system using Hamiltonian Function and Poisson Matrix. 3: 
Write the system in a cylindrical in velocity coordinate system using the formula giving how 
the Poisson Matrix and the Hamiltonian Function is transformed by change of coordinates. 
4: Make another change of coordinates in order to have the Poisson Matrix form allowing 
the application of the Key Result (Theorem 2.1). 5: Make a last change of coordinates, 
leaving the Poisson Matrix form unchanged and leading to a Hamiltonian independent of 
the last variable. 
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Panorama - As illustrated in Figure [s) the method to build the desired change of 
coordinates is made of 5 steps. The first one consists in checking that the Dynamical 
System (|9| - (10) is well Hamiltonian. This is symbolized by arrow 1 in the top picture 
of the Figure. Once this is done, we can go back into the Usual Coordinate System but 
knowing that the system writes as in the square which is in the top-left of the bottom 
picture, i.e. involving a Poisson Matrix and the gradient Vx,v^rf ^ Hamiltonian 
Function. It may be written in that form in any Coordinate System, and formula give 
how to transform the Hamiltonian Function and the Poisson Matrix while changing 
of Coordinates. The goal of the third step is to introduce a Cylindrical in Velocity 
Coordinate system which is known has being close to the Gyro-Kinetic Coordinate 
System. Writing the system in this coordinate system uses the formula giving how the 
Poisson Matrix and the Hamiltonian Function are transformed by change of coordinates. 
(At this level, it would be also possible to make a position change of coordinates fitting 
Tokamak geometry.) In the fourth step, we make another change of coordinates in order 
to have the Poisson Matrix form allowing the application of the Key Result (Theorem 
2.1 ). In the fifth step, we make a last change of coordinates, leaving the Poisson Matrix 
form unchanged and leading to a Hamiltonian independent of the last variable. 



3 Hamiltonian System 



We make here the first step: we check that Dynamical System ([9]) - ( |1Q| ) is Hamiltonian. 
Essentially, this means nothing but that there exists a coordinate system in which it 
writes: 



or 



(dq,\ 

dt 



5Vq,pi^,^ with 5 = (_° q3 



(19) 
(20) 

(21) 



where is a function. 

The Canonical Coordinates - If we use ^ and ([8| in (|9| - ([lO|, we get that in 
the Usual Coordinates (x, v) = (xi, X2, X3, '^1, '^2, '^3), trajectory : (X(t; x, v, s), V(t; x, v, s)) 
is solution to 

(22) 



dt 
dV 

~dt 



-V[$o(X)+r?$i(^)] + Vx ^^i^. 

?7 s 



(23) 



We recall that (X, V) = (Xi, X2, X3, Vi,V2, V3). 



We will show, in the next paragraph, that in coordinate system (q, p) = {qi ,q2,<l3,Pi,P2,P3) 
defined by 

A(x) 



q = x, p = v 



(24) 
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with reverse transformation given by 



X = q, V = p 



A(q) 



(25) 



the trajectory (Q(i; q, p, s), P(i; q, p, s)) ((Q,P) = (Qi, Q2, Q3, -Pi, -P2, -P3)) which is 
given by 



Q = X, P = V 



A(X) 



X = Q, V = P 



A(Q) 



is solution of a Dynamical System of the form (|T9| - (|20| or pT|) , with 

1 



A(q) 



$o(q) + ??*!(-)• 



(26) 



(27) 



Those (q, p) coordinates will be called Canonical Coordinates. 

Check of Canonical nature of Canonical Coordinates - Making calcula- 
tions suggested by (19) - (20) with given by (27) gives 



dP (VA(Q))^/ A(Q) 



-VqfL^(Q,P) 



(p-^)-V[$o(Q)+r?*i( 



dt ^ 
which applying the following formula 

(VA)^(p - A) = (VA)(p - A) + (p - A) X (V X A), 

yields 



^ P A(Q) 

dt £ ' 

dP _ (VA(Q)) A(Q) 

dt e 



A(Q)\ VxA(Q) 



V[i>o(Q) + r/$i( 



Now, using (26) we get 



dX 

dt 



dY 

'dt 



V, 



dP (VA(Q)) /dQ\ _ 
\dt) ~ 



A(Q) 



dt e 
A(Q)\ VxA(Q) 



dt 

V[*o(Q) + ??*i(-)] 

77 



(28) 
(29) 

(30) 
(31) 

H 

(32) 
(33) 

(34) 



V X Z^iAiZ^ _ v[$o(X) + r,$i(^)]. 



which is (22) - (23), and then proving the equivalence between the two systems and 
leads to the conclusion that the Dynamical System we take interest in is Hamiltonian. 
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As by-products : Poisson Matrix, Poisson Bracket, Change of Coordi- 
nates Formula - To end the first step we now list some properties we can moreover 
deduce from the fact that Dynamical System ([9| - (10) is Hamiltonian. 

In any coordinate system r = (ri, r2, rs, r4, rs, re), the Dynamical System writes: 

^ = P(R)VriI(R), (35) 

for a Poisson Matrix V, which is antisymmetric, and a Hamiltonian Function H. 

We can consider the Poisson Bracket which is defined for two regular functions / 
and ^ : ^ R by 

{f.g}{r) = (Vr/(r)) • (7^(r)(Vr^(r))). (36) 

If one of the two function is vector- valued, i.e. if for instance f : R^ ^ R^ and ^ : R^ ^ 
R, then it is defined, for any component i = 1, . . . , 6, by 

{{f,g}{r)h = (Vrf,(r)) • (P(r)(Vr5(r))). (37) 



Then, if we introduce the Coordinate Function I, defined by I(r) = r, then (35) reads 
also: 

— = {I,H}{R). (38) 

dr 



Remark 3.1. In papers of physicists, formula (38) is generally read = {r^H}. 

In another coordinate system f = (ri, f2, fa, f4, fs, fe) which is tied with the first one 
by change-of-coordinates formula f = p(r) and r = p(f) = p~^(f), the system writes 

^ = V{R)V,H{R), (39) 

where R = p(R), with Hamiltonian Function and Poisson Matrix given by 

H{i) = Hip{r)) and (P(f))i,- = {p,, pj(p(f)). (40) 



Remark 3.2. From ( |27| ) and applying ( |4Q[ )^ in Usual Coordinates, Hamiltonian Func- 
tion expression is 



leading to 



fi;,-(x,v) = i|v|' + $o(x)+r?$i(J), (41) 



and Poisson Matrix expression is 



/; 



3 



P^^(x,v) = ( (v^(x))^-(VA(x)) I ■ (43) 
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This can he easily checked calculating 



wr^ . N ^ (V^(x))^ - (V^(x)) 

7] e 



(44) 



using (29), and comparing with (22) - (|23[ ). 

Expression (41) is consistent with (40), let us check that it is the same for ( [43| ). To 
do this, let us write the change of coordinates in the following way: (x, v) = t;(q, p) 
and (q, p) = 7r(x, v). Then, for instance: 1^5 (q, p) = P2 — '^6(^7 p) = Ps — 

and 



Vi;5(q, p) 




1 

V J 



I3 
-h 



, Vt;6(q,p) 






V 1 ) 



(45) 



Hence, since {i;5,i;6} = (Vi^s) • {SVvq), expression of {V^)^q, obtained applying (40), 
is 



{t;5,t;6}(7r(x,v)) = i(^(x)-^(x)), 



(46) 



which is also the expression of {V^)^q obtained applying (43). (This may of course be 
led for the other entries of the matrix also.) 



4 Cylindrical Coordinates 



Cylindrical Coordinates in velocity - We now turn to the second step which 
consists in setting the expression of Dynamical System ([22| - ([23|) (or ([9| - ([To| or ([T^ 



(20), (27)) in a Cyhndrical in velocity Coordinate System {:>c^v\\^v±^0) which is such 
that 




Figure 6: Helicoidal trajectory and Cylindrical Coordinates for the velocity variable. 
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B 



V\\ = v ■ 



B \ B 



(47) 



This means, as shown in Figure |6] where a hehcoidal trajectory induced by a Magnetic 
Field pointing upward, that v\\ is the projection of the velocity on the direction of the 
Magnetic Field, v± is the norm of the projection of the velocity on the plan orthogonal 
to the Magnetic Field and is the angle that this projection makes in this plan. 

Hamiltonian Function and Poisson Matrix in Cylindrical Coordinates - 

Applying formula (41), we get the expression of the Hamiltonian Function in this 
System: 



Z Tj 



(48) 



We can also get the expression of the Poisson Matrix, which is heavy, applying formula 
(41). 

Here, we give its expression in the case when 



B(x) 




6 > 0, 



(49) 



only. It is: 



/ 











Kx) 

£ 


























- sin(l9) 

- cos((9) 


cos(6') 

sin^^ 

v± 




£ 











$$ 


$ 







sin((9) 


cos((9) 


-$$ 





_fc(x) 

£V± 


\ 





cos (6') 


sin(6>) 


-$ 







v± 


v± 




£V± 



(50) 



where 



^(|^(x) 

e 0x2 



db 
dx?, 



(x)) and $$ = ^ (sm(l9) -— (x) + ——- -— 

ev± 0x2 cos{0) oxs 



(x)). 



(51) 



In more general cases, we do not give the Poisson Matrix expression. Yet we mention 
a very important fact which is that 



|B(x)| 

ev± 



>0, 



(52) 



is always true. This is Important for the Darboux Method which will be reached in a 
few lines. We will denote: 



|B(x) 

ev± 



(53) 
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5 Darboux Algorithm 



Darboux Algorithm Target - The third step is the appUcation of a mathematical 
algorithm, so called the Darboux Algorithm, to build a Coordinate System (y,iX||, k^O) 
in which the Poisson Matrix has the required form, given by (15), to apply the Key 



Result. In fact, in order to manage the small parameter e, we will build the Coordinate 
System (y, u\\, /c, 0) in order to get with the following form: 



M 



0---0 

V 0---0 







(54) 



/ 



We introduce the following notations to manage change-of-coordinates mappings: (y , liy , 
k,e) = T(x, v^,6>) and {^,v\\,v±,0) = |(y, ^xy , /c, 6>), (| = T"^). 
An important and constitutive fact in the Darboux Algorithm is that the 6>- variable is 
left unchanged. 

Now, Since 

{V;{y,v\\,k,0)),j = {Ti, Y,}(C(y,«||,A:,e)),{Y„T,} = (VT,) ■ (P^(VT,)), (55) 



the bottom-right of form given in (54), results from: 

{T6,T5} = -i. 



(56) 



Remark 5.1. If we write Xq = Xq and T5 = T^^ ( [56| ) may be also read {T^*, X^} = — ^. 

In articles of physicists, this last equation reads {O^k} = — K 

In the same way, the fact that the two last lines (or columns) contain only zeros 
results from: 



{Ti,T5}=0, {Ti,T6} = 0, 

{T2,T5}=0, {T2,T6} = 0, 

{Y3,T5}=0, {Y3,Y6} = 0, 

{Y4,T5} = 0, {Y4,Y6} = 0. 

Remark 5.2. Using the same conventions as in Remark\5.1\ (IstI may also read 



(57) 



{Yj,,,Yfe} = 0, or{yi,k} = 0, {Y„,,Y4 = 0, or{yi,0} = 0, 



{Ty^,Tk} = 0, or{y2,k} = 0, 
{Yj^3,Tfc} = 0, or {ys, k} = 0, 



{Y„„Y4 = 0, or{y2,e} = 0, 
{Yy3,Y,} = 0, or{y3,e} = 0, 



{Y,,„Yfc} = 0, or{v\^,k}=0, {Y„„,Y,}=0, or {t;,,, ^} = 0. 



(58) 



Equations ( |56[ ) and ( |57[ ) are hyperbolic PDEs that need to be solve to get change- 
of-coordinates mapping T. 

First equation processing - Since (and this is a consequence of the fact that the 
6>- variable is left unchanged in the sought change-of-coordinates) 



VT6(=VT^) = (0,0,0,0,0,1)^, 



(59) 
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we deduce {T6,T5} = (VTe) • {V^{\/T^)) is the last component of {V^{\/T^)). Hence, 
equation (57) reads 



i'lT^- h r -ts— h i^ll — h UJ- 



dxi 



dxs 



dv\\ 



dv± 



(60) 



where Fn and cu are functions of (x, v\\ ,v±,0). Function uj is given by (53) and is positive. 
Solving (60) will give T^{:>c^v\\^vj_/6) which is the expression of component k of the 
Darboux Coordinates in terms of the Cylindrical in Velocity Coordinates (x, '^y , v^, 
To solve this equation we will use the Method of Characteristics. 



Method of Characteristics - Dividing by cj, which is possible since u is positive, 
(60) gives 



ST5 v^dTs v^dT^ v^dT^ v±F\\ gT5 ^ vi_ 
dv^ IBI dxi IBI dx2 IBI dxs IBI dvu IBI 



(61) 



In order to get a solution to this equation, we need to add a boundary (or initial) 
condition in a point where the involved vector field (£:v^/|B|)(Fi, F2, F3, Fy) does not 
vanish. Hence we set 



\-5\v±=u 







for a small u > 0. 



(62) 



We consider the following characteristics (Ar,V||) = (A'l, A2, A3, Vy): 



^i('^±; X, I'll , ii^) such that 
^2{v±;^^v\\^u±) such that 
^3('^±; X, I'll , ii^) such that 
Vii {v± ; X, I'll , ) such that 



SA'i ^^Fi(A'i,A'2,A'3,V||,^^,^) 
e- 



dv±_ 

dXs 

dv_i 



= £ 



\B\{X,,X2,X^) 
v^F 2{XuX2,X:„V\\,v^,e) 

IBKA*!, ^"2, A3) 
v^F3{X^,X2,X3,V\\,v^,9) 

|B|(A'i,A'2,A'3) 
v±F\\{Xi,X2,X:i,V\\,v^,e) 
\'B{Xi, X2, X3) 



X2(UI_) = X2, 
X3{u^_) = X3, 



V||(w±) = W||, 



(63) 



and the solution to (61) - (62) is given by 



nV 

^5(x, V||, v±,6>) = T5(Ar(i/;x, V||, v^), V||(i/;x, V||,vx),^,6>)+ / 

J V 



|B(Ar(5;x,^l|,^^))| 

s 



ds 



V |B(Ar(5;x, V||, v^))| 



ds. (64) 



(The last equality is gotten because of ([62|).) 

Rewriting system (63) in the following compact form, 

'AT' 



d 



Vi, 



sF(A'i,V||,^^,^), \ {u^_) 



(65) 



with F = (v^/|B|)(Fi, F^^Fs^ Fy) and applying the formula that gives the expansion of 
the solution of a dynamical system in terms of its parameter, we get that 
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x\ _ fx 

Vii j " Ivi 



^ (u_l;x,U||,w±), 



(66) 



writes 



2. ,2 



(„'||)+^F(x,«||,«^,^) + ^£f(F)(xi,«||,«^,0) 



3. ,3 



± /^2 



3! 



/:^(F)(xi,^l|,^^,^) + ... (67) 



where £f(F) is the the Lie Derivative of vector field F in the direction of vector field 
F, £|(F) = £f(£f(F)) and £|(F) = £f(£|(F)), and so on. 
Rewriting this expansion, we get 

X = ^ + £v^X^ ^e'^v'ix'^ +e^v\x'^ + ... (68) 
V|| = I'll + evi_Vl + e^uiVf + £^t;lV|f + . . . (69) 

with 

fX^\ ( X^\ 1 ( X^\ 1 

( VM =F(x,t'||,V-L,6'), ( y2 1 = ^£F(F)(x,t;||,u_L,6'), ( ^3 j = -£|(F)(x, I'll, 1-^,6'),... 

(70) 



Now, injecting (69) in expression (64) on T5, we get: 



s 

T.(x..,„.,.«) = | |B(A-(,;x.„,.,„))|* = 



where T*(l/|B(x)|) is the i^^ coefficient of the Taylor expansion of l/|B(x)|. 



Formula (71) gives the expression of T5, i.e. the expression of new variable k in 
terms of Cynlindrical Coordinates (x, , 'u^, as an expansion in e. 



Remark 5.3. In (71), it is possible to choose v as small as we want. Hence, the first 
term of the expansion is essentially the Magnetic Moment v±/\'B{-k)\. 
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On other equations and on Poisson Matrix in Darboux Coordinates 



Equation (56) was processed and gave expression of k. Equations (57) can be processed 



using similar Methods of Characteristics. A special attention needs to be given to the 
fact that Ti , . . . T4 are solutions of a PDE involving T5 and of another one involving 
Te that need to be tackled together. They will give y and u\\ as expansions in e. Once 
this done, we will have 



(72) 



or, in other words, the expression of (y,ii||,/c) in terms of {:si^v\\^v_\_^0) and as an ex- 
pansion in e. 

The other terms of the new Poisson matrix V^{y^u\\^k,0) are given computing: 
= {Ti,T2}(={T,„T,J = {yi,y2}), (P,^)i3 = {Yi,T3}(={T,„T,3} = {2/i, yg}), 



Hamiltonian Function in Darboux Coordinates - We now have to com- 
pute the Hamiltonian Function in the Darboux Coordinates. For this, we build from 
expansion (72) of T, an asymptotic expansion of ^ = T~^: 



(73) 



Then, we use the expression of the Hamiltonian Function in Cylindrical in Velocity 
Coordinates: 



1 



H^%^,v\\,v^,0) = ^(^l+^i) +$o(x)+7?$i(-). 



(74) 



and the following formula, which is gotten from (40), 



^^(y, w|| ,k,0) = H^%^{y, w|| ,k,e)) = H^%^o + + ^'^2 + + ■ • ■ ) 
^ + |B(y)|A: + $o(y) + + e^i,^(y, fc, 0) + e^H2,^{y, w,, , k, i 



(75) 



Notice that for this expansion to be valid, it is necessary to have the property we men- 
tion on page[7| r] = s^~'^ for Hi > 0. 



In expression (75), there is an important fact for the setting out of the to come Lie 



Transform based Method: the first term is independent of 0. 



6 Lie Transform based Method 

Lie Transform based Method Target - As a result of the Darboux Algorithm, 
we obtained a Poisson Matrix 7^^ (y, i^y , /c, ^) with the required form to apply the Key 
Result (Theorem 2.1), but the resulting Hamiltonian Function 



H^{y^u\\,k,e) = Ho^r]{y^u\\,k) +6:i^i,^(y,ii||,/c,6>) ^ e'^H2^rj{y^u\\,k,0) + . . . , (76) 
depends on 0. And we need to make this dependency to vanish. 
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Let us notice two important facts, which are in fact hnked: The first term in the 



asymptotic expansion (75) does not depend on 6 and the Hamiltonian Function ex- 
pressed in the Cyhndrical in Velocity Coordinates does not depend on 0. The conse- 
quence of those facts is that it is certainly possible to make the last variable to vanish 
using a mapping parametrized by e and close to identity for small e. 
Moreover, we need to build this sought mapping in such a way that it does not change 
the Poisson Matrix expression. This means that, as viewed as functions, Vr^ and 
must be the same, i.e.: 

P^^(z,w;||,j,7) =p;(z,w;||,j,7) or P^"(y, U|| , fc, ^) = P^^(y, M|| , fc, (77) 

for any z,iL'||,j,7 or y,i^||,/c,6>. Yet, among changes of variables, the symplectic ones do 
not change the Poisson Matrix expression. And among symplectic changes of variables, 
are the flows of Hamiltonian Vector Fields, which are moreover close to identity for 
small values of their parameter. 

Hence, the Lie Transform based Method consists in building a change of variables, 
parametrized by e, reading 

(y,ii||, A:,6>) ^ (z,^||, ^,7) = C{e]y,u\\,k,e), (78) 

with 

C(£; y, ^^11 , A:, 6>) = (y, u\\ , A:, 6>) + eCi (y, ^|| , ^, +^^C2(y. iiy , A:, 6>) + . . . , (79) 
which is the flow of parameter e of a Hamiltonian Vector Fields. 

The way to do - Looking for Q as the flow of a Hamiltonian Vector Field means 
looking for a Hamiltonian Function: 

G{e; y , u\\ , /c, 6>) = Go(y , u\\ , fc, 6>) + ed (y , u\\ , A:, 6>) + e^G2 (y, iiy , fc, 6>) + . . . , (80) 

such that ^(e; y, u\\ , /c, 0) is solution to: 

|^=^;VG, C(^ = 0;y,ix||,fc,^) = (y,^x||,^,^). (81) 
Hence the target becomes to find functions Go, Gi, G2, • • • , such that: 

H^{'^,w\\,j) = H^{X{z,w\\J,-f)) where A = C"^ (82) 

A result - We will use the following theorem 
Theorem 6.1. For any two Hamiltonian Functions: 

H^{y^u\\,k,e) = Ho^r^{y,u\\,k) + £:i^i,^(y, ixy , /c, 6>) ^ e'^H2^r^{y,u\\,k,0) + . . . , (83) 



H^i^,W\\J) = Ho^rj{z,W\\J) ^£Hi^rj{z,W\\J) + 7^2,77 (z, ^|| , j) + . . . , (84) 

expressed in two different variable systems (y^u\\^k^O and z^ww^j^j) , with the property 
that 

^0,77(2,^11, j) = 7^0,77(2, ^||,j), (85) 
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There exists a Hamitonian Function writing 

G{£; y, u\\ ,k,0) = Go(y, u\\ ,k,0)^ sGi (y, u\\ ,k,0) ^ e^G2{y, u\\,k,0) ^ . . . , (86) 

such that the solution C, of the Hamiltonian Dynamical System associated with G and 
with parameter €, i.e. solution to: 

^=^;VG, C(^ = 0;z,^||,j,7) = (y,ii||,fc,^), (87) 

is such that 

H^{^,w\\J) = H^{X{z,w\\J,j)) where X = C^. (88) 
Moreover, the Hamiltonian Function ^(e; z, j, 7) such that X is solution to 

= V^Vg, X{e = 0;y, u\\,k,0) = (z,^||, ^,7), (89) 

writes 

g{e]z,w\\J,j) = ^o(z,^||, j,7) + ^^i(z, J, 7) + ^^^2(z, ^y, j, 7) + • • • , (90) 
where go, gi, g2, ... are given by 

{^'0,^0,77} = Oo{Ho^r^), {gi.Ho^r]} = Oi (^0,77, ^1,77, 77, 5'o), 

{^2,^0,77} = O2(^0,77,^l,77,M,77,^2,77,^2,77,^0,5'l)7 (91) 

for differential operators Oq, Oi, O2, ... defined by recursive formula. 

Remark 6.1. Despite we do not explicitly write the dependency of X, go, gi, g2, . . 
A and Go, Gi, G2, ... with respect r] to all those functions do depend on rj. Here again, 
to be able to write the expansions with respect to e, it is necessary to have the property 
mentioned on page\^ 77 = s^~'^ for k > 0. 



The Lie Transform based Method - As in the case of the Darboux Algorithm, 
we do not do the computations that give ^ and A. We only give the steps that allows 
us to get them. The Lie Transform based Method may be summarized as: 

1. Fix: i^i,^(z,^||,j), i^2,r7(z,^||,j),•••• 
They can be fixed, a priory, with no restriction. Nevertheless those functions are 



involved in the PDEs of (91). Hence, they need to be chosen in a way that leads 



to PDEs which are as simple as possible to solve. 



2. Solve PDE in (91) to get go, ^1, ^2 recursively. 



They are quite complex but we can solve them using the Method of Characteristics. 



3. Get C solving (87). 



It is gotten as an expansion in e. 

4. Compute expression of the Lie Variables: (z, w\\ , j, 7) = (^(y, u\\ , /c, 0) 

5. Since Hi ^^^ ^2,77, • • • were fixed in the begining of the process, it is not necessary to 
compute A to get: H^{z^ w\\ , j) = i^^(A(z, w\\ , j, 7)). Nonetheless, A may be gotten 
as an expansion in e by inverting expansion of 
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7 Brief conclusion 



At the end of the day, in the Lie Coordinates System, the trajectory (Z, VKy, J, F) of a 
charged particle is solution to a Dynamical System which has the following form: 

dT 

— = Something complicated, (92) 
1=0, (93, 



dt 



Something independent of T. (94) 



This form is essentially the same as the one of (11). 



As a conclusion, we just say that the method which is summarized here is the one 



that allows us to get (11) or any other Gyro-Kinetic Approximation Model. 



Insisting one more time, the important facts in system (92 ) - (94) is that the evolution 
of the components Z and W\\ of the trajectory and of its 7-component F are uncoupled 
and that component J does not evolve and then becomes a parameter. Hence, even if 
the 7-component F is not computed, it does not preclude the computation of J, Z and 
W\\ solving a collection, parametrized by J, of four-dimensional Dynamical Systems. 
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